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We investigate the zero-temperature properties of an impurity particle interacting with a Bose- 
Einstein condensate (BEG), using a variational wavefunction that includes up to two Bogoliubov 
excitations of the BEG. This allows one to capture three-body Efimov physics, as well as to recover 
the first non-trivial terms in the weak-coupling expansion. We show that the energy and quasiparticle 
residue of the dressed impurity (polaron) are significantly lowered by three-body correlations, even 
for weak interactions where there is no Efimov trimer state in a vacuum. Eor increasing attraction 
between the impurity and the BEG, we observe a smooth crossover from atom to Efimov trimer, 
with a superposition of states near the Efimov resonance. We furthermore demonstrate that three- 
body loss does not prohibit the experimental observation of these effects. Our results thus suggest 
a route to realizing Efimov physics in a stable quantum many-body system for the first time. 


The underlying few-body dynamics in quantum many- 
body systems is fundamental to our understanding of 
nature. For instance, Fermi liquid theory is a spectacu¬ 
larly successful description of interacting many-body sys¬ 
tems in terms of the dynamics of single particles, denoted 
quasiparticles [1, 2]. Likewise, tightly bound two-body 
molecules in a molecular BEG are smoothly connected 
to Gooper pairs in a Bardeen-Gooper-Schrieffer super¬ 
fluid with decreasing attraction, as has been beautifully 
demonstrated in atomic gas experiments [3-5]. It is thus 
natural to ask whether three-body bound states lead to 
any observable effects in a many-body environment. In 
particular, Efimov famously demonstrated that there are 
an infinite number of three-body bound states (trimers) 
for identical bosons with a short range resonant interac¬ 
tion characterized by an infinite scattering length [6, 7]. 
It took more than thirty years before these Efimov states 
were first observed as a feature in the three-body loss rate 
in non-degenerate atomic gases [8-12], and more recently 
through Goulomb explosion imaging of ^He trimers [13]. 
An intriguing question is whether Efimov states can have 
observable effects on a quantum many-body system, par¬ 
ticularly outside the regime of strong three-body losses. 
So far, there have been no experiments reporting such 
effects. 

Here, we show how highly imbalanced two-component 
atomic gases can provide a different, many-body signa¬ 
ture of Efimov physics. The study of population imbal¬ 
anced gases has already yielded fundamental new insights 
concerning quasiparticles in Fermi systems (Fermi po¬ 
laron) [14-17]. With the recent identification of Feshbach 
resonances in Bose-Fermi and Bose-Bose mixtures [18- 
22], the study of quasiparticles in a BEG (Bose polaron) 
with tunable interactions is now within reach. Impor¬ 
tantly, there is no reason to expect three-body correla¬ 
tions to be small for the Bose polaron, in contrast to 
the Fermi polaron, where these are suppressed by Fermi 
statistics in the two-component system [23, 24]. To have 


enhanced three-body correlations in Fermi systems, one 
requires unequal masses [25], low dimensions [26-28], or 
at least three different species [29, 30]. 

The large majority of theoretical studies of the Bose 
polaron have been based on mean-field theory [31-35], or 
the Frohlich model [36-42]. These approaches, however, 
do not include three-body Efimov correlations, and the 
same applies to recent variational [43] and field theoret¬ 
ical [44] studies. Only recently did a perturbative calcu¬ 
lation demonstrate that three-body correlations indeed 
show up at third order in the impurity-boson scattering 
length a, leading to important effects such as a logarith¬ 
mic contribution to the polaron energy [45]. 

In this Letter, we develop a variational ansatz for the 
Bose polaron that includes three-body Efimov physics, 
and in addition recovers all contributions to the energy in 
the weak-coupling limit up to and including the logarith¬ 
mic term at order . We show that Efimov correlations 
significantly lower the energy and quasiparticle residue, 
even for scattering lengths substantially smaller in mag¬ 
nitude than a_ < 0, the scattering length where the deep¬ 
est Efimov trimer crosses the three-body continuum in a 
vacuum. As a approaches a_, the energy exhibits a clear 
avoided crossing, reflecting a smooth change in the quasi¬ 
particle character, from single-particle to three-body. We 
furthermore demonstrate that these effects can be ob¬ 
served before encountering resonantly enhanced three- 
body loss, which has so far been the main experimental 
signature of Efimov physics. 

Model.- We consider an impurity particle immersed 
in a weakly interacting BEG with 1, where no 

denotes the condensate density and ub the boson-boson 
scattering length. The addition of an impurity atom to 
the condensate introduces two additional length scales: 
The impurity-boson scattering length a, and the so-called 
three-body parameter. The latter sets the scale of one 
Efimov trimer in the universal limit of 1/a ^ 0, from 
which all other trimer energies can be determined. In the 
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cold atoms context, it is natural to use a_ as the three- 
body parameter. It has recently emerged in studies of 
identical bosons [46-48] that this parameter is universally 
related to the van der Waals range. Thus, we will fix a_ 
through the two-body physics. 

For this purpose, it is convenient to use the two- 
channel Hamiltonian (setting h and the volume to 1): 

^ [-Bk/3k^k + ek4ck + (ck + I'o) dl.dk 

k 

+ X] + h.c?j + 5 T] {dqCd-kbk + h.c.) , (1) 

k k,q 

which measures the energy with respect to that of the 
BEC. Here, and cj^ create a boson and the impurity, 
respectively, with momentum k and single-particle en¬ 
ergy Ck = 12m. The impurity interacts with a boson 

by forming a closed channel molecule, created by with 
dispersion = Ck/2 and bare detuning z/q. The coupling 
strength for this process, g, is taken constant up to a 
momentum cut-off A. We follow the usual weak-coupling 
prescription, relating the boson operator b to the Bo- 
goliubov excitation /3 via 6k = — '^k/^Lk with real 

and positive coherence factors = [1 + (ck + /i)/Ek]/2 
and = [—1 + (ck + /i)/Ek]/2, dispersion E\^ = 

yek(ek + 2/i), and chemical potential fi = A'Ka^no/m. 
The impurity is assumed to have the same mass m as 
the bosons, as we do not expect our results to change 
qualitatively in the mass imbalanced system. 

By considering impurity-boson scattering, we relate 
the bare molecule detuning z/q, coupling and momen¬ 
tum cut-off to the scattering length a and the effective 
range ro [49, 50]. The additional length scale ro provided 
by the closed channel fixes the three-body parameter of 
Efimov physics. In particular, by solving the impurity- 
boson-boson problem in vacuum with ub = 0, we find 
a_ ~ 2467ro, which is consistent with calculations in¬ 
volving realistic interatomic potentials for a range of het- 
eronuclear systems [51]. This large separation of scales 
means we can always take |a/ro| ^ 1, where the effect of 
the closed channel on two-body processes is negligible. 

Variational approach.- We consider the wavefunction: 

It/,) = [ aocl + akcLk/3k + \Y1 “kik^cLk^-k/kyk^ 
V k kik2 

+ 7o 4 + E 7k<iLk/3k j 1^), (2) 

with |<I>) the wavefunction of the interacting BEC. The 
first line of Eq. (2) describes the bare impurity and its 
dressing by one or two Bogoliubov modes. The second 
line describes the impurity bound to a boson and form¬ 
ing a closed channel molecule, which can be dressed by 
a Bogoliubov mode. The variational approach consists 
in taking the stationary condition {d'lpl {H — E) lip) =0 


where the derivative is with respect to any of the ex¬ 
pansion parameters. This procedure yields five coupled 
equations for the energy E [50]. Defining 

_ 9 ^ Zlk '^kO^k ^ _ 9 ^ Z^k^ '^^k^Q^kk^ 

9 ZT* 7 / ’ 77' d IT ’’ 

E -Uq E - Z/q - 6^ - Ek 

which remain finite in the limit z/q ^ oo, and carrying 
out the renormalization procedure, we finally arrive at 
the coupled integral equations 
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where = E — E\^^ —E \^2 ~^ki+k 2 - Here, the medium 

T matrix is 

T-\k„E) = ^ - ^{E -eij- n(ki, E), (5) 

with n(ki,i;) = EkJ^L/(-®“-^k2 -eki+k2)+2mr/q] 

the pair propagator in the BEC and = m/2 the re¬ 
duced mass. 

A few comments on our variational approach are in 
order. Eirst, had the two equations (4) not been coupled, 
the first of these would simply reproduce the mean-field 
polaron energy, E = Annoa/m (for a non-zero ro there 
is a negligible correction ^ O{n‘^a^ro) as ro is by far the 
smallest length scale). Second, in the limit no ^ 0 the 
second equation exactly yields the Efimov spectrum for 
the impurity and two identical, non-interacting bosons. 
This explicitly demonstrates that our variational ansatz 
contains the relevant three-body Efimov correlations. Eor 
a discussion of the three-body problem with an effective 
range, we refer the reader to Refs. [52-54]. 

Eurthermore, since our wavefunction includes all two- 
body excitations, it recovers all second order perturba¬ 
tion terms in a. This is because all relevant Eeynman 
diagrams contain at most two concurrent boson propaga¬ 
tors (neglecting noa^ corrections) [45]. Only a subset of 
these diagrams are produced if one restricts the approach 
to one Bogoliubov mode, as in Ref. [43]. Although our 
ansatz does not reproduce all third order self energy dia¬ 
grams, it does indeed capture the logarithmic term in the 
energy, since this originates from diagrams containing at 
most two simultaneous Bogoliubov excitations [45]. 

Einally, we remark that had we used a single-channel 
model from the outset, we would have arrived at equa¬ 
tions identical to Eqs. (4) with ro = 0 in Eq. (5). How¬ 
ever, in this case we would need to set the three-body 
parameter in some other way, for instance through the 
use of a Gaussian cut-off. 
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Figure 1. (a) Spectrum and (b) residue of the Bose polaron 
for = — 1 and a^ja- — —1/50. We show the results 

from our variational calculation (purple, solid), and the vari¬ 
ational wavefunction excluding three-body correlations and 
effective range effects (blue, dashed). The gray lines are the 
results of perturbation theory: Mean-field (dashed), second 
order [38] (dot-dashed), and third order [45] (dotted). We also 
display the energy of the deepest Ehmov state in a vacuum 
(red, dashed), and for a > 0 the threshold of the atom-dimer 
conti nuum giv en by the impurity-boson bound state energy 
— (y^l — 2ro/a — l)^/(2mrro) (black, solid). The insets illus¬ 
trate the crossover in the wavefunction (the atoms are not 
drawn to scale). 


Results.- The Bose polaron is characterized by its 
quasiparticle properties, i.e. the energy E and residue 
Z = |(aop [50]. In Fig. 1, we display these quanti- 
ties calculated from Eq. (4) as a function of l/{nQ a) 

for realistic experimental parameters. For weak interac- 

- 1/0 

tion, a) — 1 , the polaron energy is close to the 

mean-field value Airano/m and Z 1. With increasing 
coupling, E starts to deviate from the mean-field result, 
whereas there is better agreement with second [38] and 
third [45] order perturbation theory. This is as expected 
since, as argued above, our variational approach contains 
all diagrams contributing to the energy at these orders 
of perturbation theory. While Z = 1 within mean-field 
theory, and the second order perturbative expression [45] 
is negative in the range plotted, the agreement at third 
order with our variational ansatz is very good, leading 
us to conclude that although the variational ansatz does 
not contain all third order terms for the residue, it likely 
includes the most important ones. 


I 






Figure 2. Residue of the Bose polaron as a function of 
l/{nl^^a) and l/(nj^^|a-|) for a^ja- — —1/50. The white 
dashed line is a = a-. For sufficiently small 1/(72^^|a-1), the 
residue remains close to 1 around unitarity, implying that the 
Efimov trimer has disappeared. 


Closer to resonance, three-body correlations become 
important, leading to an avoided crossing between a state 
best described as a single particle interacting attractively 
with the medium and an Efimov trimer dressed by inter¬ 
actions, as seen in Fig. 1 and illustrated in the insets. The 
residue correspondingly decreases as the wavefunction re¬ 
sembles the three-body Efimov state. This crossover is 
of course completely missed by the perturbative expres¬ 
sions. It is also missed by the variational ansatz used 
in Ref. [43], which includes only two-body correlations, 
i.e., only the first two terms in Eq. (2). Within this ap¬ 
proach the quasiparticle residue only goes to zero once 
the polaron approaches the atom-dimer continuum. Our 
analysis is the first to include Efimov physics for the Bose 
polaron, which we see strongly changes its properties at 
a qualitative level. Even outside the regime of Efimov 
physics in vacuum, there are strong effects of Efimov 
physics: For instance, at a = a_/2 the magnitude of 
the energy is increased by 50% when comparing with the 
variational approach restricted to two-body correlations. 

To investigate the strength of the coupling between 
the two branches involved in the avoided crossing, we 
show in Fig. 2 the residue as a function of l/(ny^a) 
and 1/(71^^Ia_|). Focussing first on the region where 
l/(nQ^ l<^-l) ^ I 5 we see that the crossover region is lo¬ 
cated around a ^ a_, where Z decreases to almost zero. 
Furthermore, the range of scattering lengths over which 
this transition occurs decreases with decreasing |a_|, re¬ 
flecting how the coupling to the Efimov state becomes 
weaker. This can be understood from a simple dimen¬ 
sional analysis as follows: Treating momentarily the Efi¬ 
mov trimer as a point particle, the coupling of this state 
to two bosons and the impurity involves four particles, 
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three in-going and one out-going or vice versa. The cou¬ 
pling strength for this process must therefore scale as 
Length/Mass. In the avoided crossing region, there is 
only one relevant length scale, a ^ a_, and consequently 
the coupling strength ^ a_/m. On the other hand, in the 
regime l/{nQ ^ 1 the perturbative energy far ex¬ 

ceeds that of the Efimov trimer in vacuum when a ^ a_ 
and thus the quasiparticle remains well-defined even be¬ 
yond this point, as seen in the figure. Eventually, of 
course, the residue is suppressed as the polaron becomes 
dimer-like for sufficiently small and positive a. 

To investigate the dependence on ub, we plot in Eig. 
3(a-b) the energy and residue of the Bose polaron for 
Uq a_ = —1/4 and three different values of the boson- 
boson scattering length. In both quantities we clearly 
observe how the crossover becomes sharper with increas¬ 
ing ae/ltt-l, since boson-boson repulsion suppresses the 
formation of the Efimov trimer. An even larger suppres¬ 
sion is expected to occur for any four-body Efimov state 
[55] consisting of the impurity and three bosons, if such 
states exist (these have so far only been investigated for 
unequal masses [56]). Indeed, the impact of four-body 
states on the many-body system will be further reduced 
since one expects them to appear for \a\ < |a_|, resulting 
in a smaller coupling to the bound state similar to that 
shown in Eig. 2. It is thus likely that such states can 
be ignored even in a regime where the coupling to the 
three-body state is strong. 

While the energy and residue of the Eermi polaron 
have been measured using radio-frequency (RE) spec¬ 
troscopy [14-16], the presence of the Efimov trimer in 
a cold atoms context has hitherto only been observed in¬ 
directly, as either a significant increase (suppression) in 
the three-body loss rate for negative (positive) scattering 
lengths [8-11], or via enhanced loss from the RE associa¬ 
tion of trimers close to the atom-dimer threshold [12]. An 
important question is therefore whether the shifts in the 
polaron energy and residue due to Efimov correlations 
can be observed before encountering prohibitively fast 
losses. To investigate this, we note that the main loss 
process - three-body recombination to a deeply bound 
molecule and an atom - is local, since it requires three 
atoms to be separated within a distance comparable to 
the size of the outgoing states. Hence, the rate of this 
loss is captured, within a perturbative approach, by the 
expectation value [57, 58] 

r = A(^| ^ 4_k,4,^'k.dp-k.lV’). (6) 

kik2P 

Making the replacement Pk at short range, we get 

T l7/cP- Note that the constant A should be 

fixed by comparison with experimental data outside the 
regime of resonantly enhanced losses; thus, in Eig. 3(c) 
we plot T/A. We see that the loss rate increases with 
decreasing Z due to the onset of Efimov correlations, and 



Figure 3. (a) Energy, (b) residue, and (c) three-body loss 

rate for n^^la-l = 0.25 and different values of aB/|a-|: 0.01 
(purple, solid), 0.1 (purple, dashed), 0.5 (purple, dotted). In 
(a) we also display the energy of the deepest Efimov state 
in vacuum (red, dashed). The blue lines are the variational 
results excluding three-body correlations as well as effective 
range effects. 


saturates to that of the vacuum Efimov state. Eigure 3, 
however, demonstrates a crucial point: The reduction in 
the energy and residue due to three-body correlations is 
still substantial away from a = a_, where the decay rate 
can be more than two orders of magnitude smaller than 
that of the vacuum Efimov state. Erom this we conclude 
that the lifetime of the impurity is indeed long enough 
in this regime to allow a detection of three-body effects. 
This mirrors previous work on trimers in a quasi-two- 
dimensional geometry, which argues that Efimov physics 
can be probed by coupling an Efimov trimer to a trimer 
extended within the plane [59]. 

Conelusions.- Using a variational ansatz which in¬ 
cludes three-body Efimov correlations as well as the low¬ 
est non-trivial terms in the weak-coupling expansion, we 
demonstrated that Efimov physics dramatically changes 
the properties of the Bose polaron. The energy and quasi¬ 
particle residue are lowered significantly, even when there 
is no Efimov trimer state in a vacuum. As the critical 
scattering length for the emergence of the trimer state is 
approached, the energy exhibits an avoided crossing and 
the residue quickly decreases, reflecting the fact that the 
impurity smoothly becomes bound in an Efimov trimer. 
Crucially, we showed that three-body loss associated with 
the Efimov state does not prohibit the observation of 
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these effects. In conclusion, our results demonstrate how 
Efimov physics qualitatively changes the properties of a 
quantum many-body system, and how it can be observed 
in signals other than loss. 
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IMPURITY-BOSON SCATTERING 


The bare molecule energy z/q is not a physical observable since the molecule is dressed by open channel impurity- 
boson pairs. Likewise, the coupling is assumed to be a constant g for momenta below an ultraviolet cut-off A, which is 
not related directly to any observable. We can relate these parameters to observables by considering impurity-boson 
scattering via the closed channel molecule. In a vacuum, the scattering matrix is [49, 58] 


"7vac(^) — 


^ ^0 5'^nvac(0) 5^^^nvac(^) 


(SI) 


where nvac(<^) = “ k‘^/2m is the vacuum pair propagator with energy cj, = m/2 is the reduced 

mass, and Anvac(<^) = nvac(<^) ~ nvac(O)- Equating (SI) evaluated at the scattering energy uj = with the 

effective range expansion 7^ac(p^/2mr) = 27Ta/mr{l — k?are^/2 + iko), where ro is the effective range, we obtain 
a = — z/o)/27r and ro = —27r/m^^^. In this way we replace the bare parameters z/q and A with the 

low energy observables a and ro for the impurity-boson scattering. Then we are at liberty to take the limits z/o ^ oo, 
A oo. 


VARIATIONAL EQUATIONS 

Taking the derivative {K — E) \^Ij) = 0 with respect to Ofk, Q^kiks? To, and 7k yields the five equations 

Ea^ ~ 3^ ^k7k, 

k 

(E - ek - Ek) Qfk =9U\^Uo + ^v^^Tk, 

-E'kik2^kik2 ~9 ('Tki'^k2 E 7k2'^ki) , 

{E - z/o)7o =9V^^o E9^ ^k<ak, 

k 

[E -e^-z/o-Ek) 7k =^V^<^k -^'^k<ao + ^ ^ ^^k'Q^kk', (S2) 

k' 

where = E - E]^^ - E]^^ - eki+k 2 - 


QUASIPARTICLE RESIDUE 

Assuming that the variational wavefunction in Eq. (2) is normalized, = 1, the residue is given by Z = |(aoP- 

Carrying out the renormalization procedure, we arrive at 

f Vmr] 

_ V ^ ^k ^ ; _ 

f Vn^V _ A ^ _ 1 _ Zp J_ ( V^^k + ^kr ^ ^ | ik _J_ (^kiCk2+^k2Cki ) 

E 2^k E j ' ^ A-^k E—e^^ — E\^ j A-^k g‘^ Z^kik2 


(S3) 



DECAY RATE 


In the decay rate, we need to evaluate |7kP- This is evaluated similarly to Z above: 

V ^ 

I |2 -^k ^2 

, f y/^'n _ A ^ I I ^ ^ ^ y^ (^kiCk 2 +^k 2 Cki) 

£; Z^k £; j ' 5-2 “T ^ £;_£im_^^ j n“ ^2 “r Z^k^k2 



